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On the Constancy of the Photon Index of X-ray spectra of 
4U 1728-34 through all spectral states 

Elena Seifina 1 and Lev Titarchuk 2 
ABSTRACT 

We present an analysis of the spectral properties observed in X-rays from 
Neutron Star X-ray binary 4U 1728-34 during transitions between the low and 
the high luminosity states when electron temperature kT e of the Compton cloud 
monotonically decreases from 15 to 2.5 keV. We analyze the transition episodes 
from this source observed with BeppoSAX and Rossi X-ray Timing Explorer 
(RXTE) satellites. We find that the X-ray broad-band energy spectra of 4U 1728- 
34 during all spectral states can be modeled by a combination of a thermal 
(black body-like) component, a Comptonized component (which we herein de- 
note COMPTB) and a Gaussian component. Spectral analysis using this model 
evidences that the photon power-law index Y is almost constant (r = 1.99 ±0.02) 
when kT e changes from 15 to 2.5 keV during these spectral transitions. We ex- 
plain this quasi-stability of the index V by the model in which the spectrum is 
dominated by the strong thermal Comptonized component formed in the transi- 
tion layer (TL) between the accretion disk and neutron star surface. The index 
quasi-stability takes place when the energy release in the TL is much higher than 
the flux coming to the TL from the accretion disk. Moreover, this index stabil- 
ity effect now established for 4U 1728-34 during spectral evolution of the source 
was previously suggested for a number of other neutron binaries (see Farinelli & 
Titarchuk, 2011). This intrinsic property of neutron star is fundamentally dif- 
ferent from that in black hole binary sources for which the index monotonically 
increases during spectral transition from the low state to high state and saturates 
at high values of mass accretion rate. 
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Introduction 



The evolution of spectral parameters of compact objects in X-ray binaries is of great 
interest for understanding the nature of compact objects. It is well known that a number 
of black hole (BH) candidate sources demonstrate correlations between their 1-10 Hz quasi- 
periodic oscillation frequencies (QPOs) v L and photon power- law index T during spectral 



transi tion when sources evolve from the low to high states [see IShaposhnikov fc Titarchuk 



( 120091 ). hereafter ST09]. Then the definition of the spectral state is related to the level of 
soft blackbody emission presumably related to mass accretion rate. In the high states of BHs 
these index-QPO frequency correlation s sometimes show a saturati on of T at high values of 
u L . On the other hand ST09 [see also iTitarchuk fc Seifina I (120091 ). hereafter TS09] found 
that T saturates with mass accretion rate in almost any case of a BH binary. This saturation 
effect can be considered as a black hole signature or signature of the converging flow into 
BH (ST09 and TS09). The question naturally arises how the spectral index behaves as a 
function of mass accretion rate or as a function of cutoff energy of the spectrum in neutron 
star (NS) sources. 



Recently iFarinelli fc Titarchuk I ( 1201 ll ). hereafter FT11, collected X-ray spectra ob- 
tained by BeppoSAX from quite a few NS sources Sco X-l, GX 17+2, Cyg X-2, GX 340+0, 
GX 3+1 and GS 1826-238. Their results probably indicate that the value of the photon 
index slightly varies around 2 independently of the spectral state (or electron temperature o f 
Compton cloud ) at lea st for this particular sample of NS spectra [see lDi Salvo et al. I ( 120001 ) : 
Farinelli et al. I ( 120081 )]. However, the available data for those sources were taken when these 
sources were in the high state or in the low state but nobody has analyzed up to now the 
spectral evolution from the low state to the high state for some particular NS source. 

A suitable candidate for the study of the spectral evolution in NS is the so called atoll 
4U 1728-34 which exhibits a remarkable spectral transition from the low state to the high 
state and vice versa. 4U 1728-3 4 (GX 354—0) was fi rst r esolved by UHURU scans of the 



Galactic center region in 1976 [see lForman et al. I ( 119761 ) and lBradt et al. I ( 11993)1 . Then type 



I X-ra y bu rsts from 4U 172 8 -34 w ere discovered during SAS-3 observations by iLewin et al. 
(119761 ) and lHoffman et al. I ( 119761 ). Further the bursting behavior was subsequently studied 
in detail u sing extensive ob s ervat ions by SAS-3, which accumulated 96 bursts in total. Using 



these data lBasinska et al. I (11984 ) presented evidence for a narrow distribution of peak burst 
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fluxes, as well as a correlation between the peak flux and the b urst phase. The dist ance 
to the source in the range of 4.2 - 6.4 k pc has been estimated b y Ivan Paradijs I (119781 ) and 
confirmed by lBasinska et al. I (119841 ) and lKaminker et al. I ( 119891 ) using measurements of the 
peak burst fluxes. 

A radio counterpart of 4U 1728-34 was detected durin g VLA (at 4.86 G Hz) observations 
with a variable flux density in the range of ~0.3 - 0.6 mJy (IMarti et al.lll998l ). The estimated 
extinction of the source is Ay ~14 and a precise position following from the det ection of the 
radio counterpart allowed to identify this source as a K = 15 infrared source ( IMarti et al. 
19981 ). Long-term Ariel-5 measurements, as well as the extensive monitoring by the All- 
Sky Monitor (ASM) on th e board of RXTE , suggest the presence of a long-term quasi- 



periodicity about 63-72 d ( jKong et al. I Il998f ). RXTE/PCA observations of the source in 



1996 led to the discovery of nearly coherent millisecond oscillations during the X-ray bursts 



(jStrohmayer et al. Ill996h . 



Titarchuk fc Osherovich I (119991 ) presented a model for the radial oscillations and diffu- 



sion of the perturbation in the transition layer (TL) surrounding the neutron star. Using 
dimensional analysis, they have identified the corresponding radial oscillation and diffusion 
frequencies in the TL with the low-Lorentzian ul and break frequencies V\, for 4U 1728-34. 
They predicted values for i/& related to the diffusion in the boundary layer, that are consistent 
with the observed u^. 

Subsets of the bur sts ob served during the PCA observations have been also studied 
by Ivan Straaten I (120011 ) and iFranco I (120011 ) with particular attention to the relationship 
between the appearance of burst oscillations and a value of mass accretion rate. 



Titarchuk fc Shaposhnikov I (120051 ). hereafter TS05, analyzed RXTE/PCA observations 



of 4U 1728-34 in energy range from 3 keV to 40 keV and they found that using the model 
comprising two Comptonization components (BMC0) the photon index T is consistent with 
being quasi-constant (around 2.2), while the low quasi-periodic oscillation (QPO) frequency 
does not exceed 10 Hz, then T monotonically incr eases to values of 6. Mo reove r, using broad- 
band observations of 4U 1728-34 by BeppoSAX. |Pi Salvo et al. I (120001 ) and iPiraino et al. 
(120001 ) fitted the X-ray spectra of 4U 1728-34 in a wide energy range from 0.1 keV to 200 
keV by a sum of Blackbody component pl us a therma l Com p tonization spectrum, usua lly de- 
scribed by the XSPEC COMPTT model jTitarchuk I Jl994f ). lHua fc Titarchuk! (jl995h ]. The 
TS05 model cannot fit the BeppoSAX data and moreover this model has more parameters 



X BMC is so called "Bulk Motion Comptonization" XSPEC Model (see detailes in 
tion 6.2.10 of "User's Guide of an X-Ray Spectral Fitting Package XSPEC v.12.6.0" 



sec- 
and 



http : / /heasarc . gsfc . nasa. gov /xanadu /xspec /manual / Additive . html ) 
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than that of model of iDi Salvo et al. I (120001 ). IPiraino et al. I fl2000f ) and FT11. Naturally 
we pose the question if the FT11 model can fit the BeppoSAX along with RXTE data and 
what kind of the dependence of index vs spectral state can be found. 

In this Paper we do try to answer a fundamental question on the possibility to distinguish 
BH from NS systems, which extensively attempted to be solved without considering the mass 
of the compact object as a main argument, in particular without using optical counterpart 
data to measure mass function. In this way some methods have been proposed to identify 
systems contains BHs using X-ray observational properties only. The strong rapid variability 
was firstly considered as a signature of the presence of BH (Oda et al., 1971), until the 
same rapid variability was detected in accreting NSs (ITennant et al. Ill986l ). Now it is well 
established that Galactic BH candidates (BHCs) demonstrate tw o spectral states, the "high 



state " (HS) and "low state" (LS) and transition between them (iRemillard fe MacClintock 



2006 ). Howeve r, sometimes so cal led atoll-NS sources § also show the "high" and "low" 
spectral states |D' Ai et al. I (120061 ) . TS05]. Therefore this property requires more detailed 
investigations BHs vs NSs. Specifically, the HS spectra of BHCs are characterized by thermal 
emission at ~1 keV presumably originated in the accretion disk, along with a steep power 
law tail whose photon index (T=2— 3) monotonically increases with mass accretion rate (see 
ST09). While the LS spectra show much weaker disk emission than that in the HS spectra 
and harder power-law tail (which photon index is around T ~ 1.7). This hard component 
generally believed is a result of thermal Comptonization of soft (disk) photons in a hot gas 
(Compton cloud) in vicinity of the compact object. 

It is worth to point out that BHs, in contrast with NSs, sometimes demonstrate more 
complicate X-ray spectrum. For example, the RXTE spec tra of BH GRS 1915+105 re quire 
two Comptonization components, soft and hard ones [see iTitarchuk fe Seifina I (120091 )]. In 
this case one can clear see the evolution of two p hoton indices = 1.7—3.0 a nd T 2 = 2.7—4.2 



for the hard and soft components respectively. ITitarchuk fe Seifina I ( 120091 ) argued that the 
index saturation effect of the hard component is due to Comptonization of the soft (disk) 
photons in the converging flow into BH and that of the soft component is due to the thermal 
Comptonization in the transition layer when mass acc retion rate increases. These c onclusions 
were later supported by Monte-Carlo simulations by lLaurent fe Titarchuk I (1201 if ). 



Moreover, in the description of BH and NS X-ray low state spectra with the thermal 
Comptonization model, there is essential difference between these types of the compact 
sources. The electron temperature of the Compton (scattering) cloud kT e is usually lower 



2 Here, we use a term of the atoll-NS sources to specify NS X-ray binaries characterized by a specific 
'e" -shaped track in color-color diagram. 
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for NSs, kT e < 25 keV than that for BHs, kT e > 50 keV [see Ichurazov et al. I (119971 )]. The 
lower electron temperature in NSs than that in BHs is a consequence of the additional cooling 
provided by the NS surface whi ch reflects X-ray photon s and ultima tely determines a value 
of the CC electron temp erature Titarchuk et al. I ( 119981 ) and see also ISunyaev fc Titarchuk 
jl989( l. lKluzniakNl993h ]. 

This fact of t he observed difference betw een the black hole and neutron stars was recently 
also discussed by iReynolds fc Miller I (120111 ) . They conclude the observable difference of kT e 
in these types of the sources is as an evidence of the absence and presence of a solid surface 
in BHs and NSs respectively and this fact can be considered as an indirect evidence f or the 
exis tence of the event horizon in BHs. It is worth to point out that iTitarchuk et al. I (119981 ) 



and ITitarchuk &: Fiorito I (120041 ) previously came to the similar conclusions based on the 
analysis of the Compton cooling of X-ray emission region (transition layer) in the presence 
(NS) and absence (BH) of the reflection surface. 

Thus, the basic property, which differs BHs from NSs, is the presence of event horizon as 
well as a converging flow in vicinity of BH (jEbisawa et al. 1119961 ). In fact, close to the event 
horizon, the strong gravitational force dominates the pressure forces which leads to an almost 
free fall converging flow of accreting material into a BH. The dynamical Comptonization of 
low energy photons off fast-moving electrons dominates the thermal Comptonization at high 
mass accretion rate, then the plasma temperature of converging flow is less that 10—1 5 



keV, and as a result an exte nded steep power law is formed [ see ITitarchuk et al. I ( 119971 ) 



Titarchuk fc Zannias I (119981 ) and lLaurent fc Titarchuk I (119991 )]. This kind of the spectra 



are observed in the soft state of BH binaries (see e.g. ST09 and TS09). 

On the other hand in NS sources the radiation pressure forces become dominant close to 
their surface and thus a free fall should be suppressed at high mass accretion rates. Does the 
presence of the firm surface in NS makes any difference as for the dependence of the photon 
index T vs mass accretion rate with respect to that established in BHs (see ST09, TS09 and 
ST10)? Furthermore can the index saturation detected in many BHs with mass accretion 
rate exist only i n BH sources and cannot be observed so far in NS sources? For example 



Di Salvo et al. I ( 120061 ). studying of low mass X-ray binaries hosting neutron stars, concluded 



that it is unlikely to distinguish BHs from NSs based on their X-ray spectra. However FT11 
argue that in NS sources the index T weakly varies until the soft photon illumination of the 
transition layer (TL) Qd is much smaller than the energy release in the TL, Q cor . We try 
further to test this kind of the index behavior in the NS source using X-ray observations of 
the atoll source 4U 1728-34 and compare it, if it's possible, with the index dependence on 
mass accretion rate established in BHs. 



In this Paper we present the analysis of the BeppoSAX available observations during 1998 
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- 1999 years and RXTE/PCA /HEXTE observations during 1996 - 2000 years for 4U 1728- 
34. In §2 we present the list of observations used in our data analysis while in §3 we provide 
the details of X-ray spectral analysis. We analyze an evolution of X-ray spectral and timing 
properties during the state transition in §4 and §5. We discuss our results and make our 
conclusions in §6 and §7. 



Data Selection 



Broad band energy spectra of the source were obtained combining data from three 
BeppoS AX Narrow Field Ins truments (NFIs): the Low Energy Concentrator Spectrometer 
[LECS; iParmar et al. 1 dl997hl for 0.3 - 4 keV, the Medium Energy Concentrator Spectrome- 



ter [MECS; 



Boella et al. 



Frontera et al. 



(119971 )] for 1.8-10 keV and the Phoswich Detection System [PDS; 
( 1l997l )] for 15— 60 keV. The SAXDAS data analysis package is used for 
processing data. For each of the instruments we perform the spectral analysis in the energy 
range for which response matrix is well determined. The LECS data have been re-normalized 
based on MECS. Relative normalization of the NFIs were treated as free parameters in model 
fitting, except for the MECS normalization that was fixed at a value of 1. We checked after 
that the fitting procedure if these normalizations were in a standard range for each instru- 
ments 3 ! 



Specifically, LECS/MECS re-normalization ratio is 0.92 and PDS/MECS re-normalization 
ratio is 0.97. In addition, spectra are rebinned accordingly to energy resolution of the instru- 
ments in order to obtain significant data points. We rebinned the LECS spectra applying a 
rebinning template for grouping (lecs_2. grouping) with an energy depending binning factor 
used in GRPPHA of XSPEC0. Also we rebinned the PDS spectra with linear binning factor 
2, grouping two bins together (resulting bin width is 1 keV). In Table 1 we list the BeppoSAX 
observations used in our analysis. 



We also analyzed the available data obtained with RXTE ( IBradt et al. Ill993l ) which 
have been found in the time period from February 1996 to July 2000 [see also a review by 
Galloway et al. I (120081 )]. In our investigation we selected 127 observations made at different 
count rates (luminosity states) with a good coverage of rise-decay flare track. We have made 
an analysis of RXTE observations of 4U 1728-34 during four years for 8 intervals indicated 
by blue rectangles in Figured] (top). We have also analyzed two BeppoSAX observations 



3 http://heasarc. nasa.gov/docs/sax/abc/saxabc/saxabc. html 
4 http://heasarc. gsfc.nasa.gov/FTP/sax/cal/responses/grouping 
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which dates marked by green triangles there. 

Standard tasks of the HEASOFT/FTOOLS 5.3 software package were utilized for data 
processing. For spectral analysis we used PCA Standard 2 mode data, collected in the 3 - 
20 keV energy range. The standard dead time correction procedure has been applied to the 
data. The average dead time correction is in the range 3 - 10 % depending on the count 
rate value. 

HEXTE data have been also used in order to construct broad-band spectra. We sub- 
tracted the background corrected in off-source observations. To exclude the channels with 
largest uncertainties we used only data in 20 - 60 keV energy range for the spectral analysis. 
The HEXTE data have been re-normalized based on the PCA. Typical PCA/HEXTE renor- 
malization factor is 0.98. We used the data which are available through the GSFC public 
archive (http://heasarc.gsfc.nasa.gov). In Table 2 we list the groups of RXTE observations 
which cover the source evolution from quiescent to flare events. 

Note that we did not use any normalization factor to normalize between BeppoSAX 
and RXTE data. We also used public 4U 1728-34 data from the All-Sky Monitor (ASM) 
on-board RXTE. We retrieved the ASM light curves (in 2—12 keV energy range) from the 
public RXTE /ASM archive in MIT@. In the bottom panel of Figure [1] we show a mean count 
rate (blue dashed line) during 1996—2010 interval of ASM/RXTE monitoring observations of 
4U 1728-34. In this panel one can also see a long-term quasi-periodic variability of mean soft 
flux during ~ six years cycle. We investigate available periods of slow variability (indicated 
by green during which we also have the BeppoSAX observations of 4U 1728-34 (see also the 
upper panel of Figure [1]). 

We use definitions of the low and high luminosity states to relate these states to the 
source luminosity and we demonstrate that during the high-low state transition the electron 
temperature of Compton cloud changes from 2.5 keV to 15 keV respectively and vice versa. 
Thus the "high spectral state" corresponds to "low electron temperature state" and vice 
versa the "low spectral state" corresponds to "high electron temperature state" . During the 
flare, seen in the ASM light curve, the electron temperature kT e usually decreases from 15 
keV to 2.5 keV. We introduce a definition of a "burst" to point out a significant increase of 
X-ray flux (about factor of 5) with respect to the persistent emission level. Specifically, we 
call the "burst" when ASM count rate is greater than 10 c/s. We associate the count rate 
increase with the increase of mass accretion rate. 



5 http: / /xte. mit.edu/ ASM_lc.html 



-8 - 



3. Spectral Analysis 



In our spectral data analysis we use a model which consists a sum of Comptonization 



(COM PTB) component, [COMPTB is the XSPEC Contributed modefl, see iFarinelli et al. 



( 120081 ). hereafter F08], a soft blackbody component of temperature T BB and Gaussian line 
component. The COMPTB spectral component has the following parameters: temperature 
of the seed photons T s , energy index of the Comptonization spectrum a (= T — 1), electron 
temperature T e , Comptonization fraction / [/ = A/(l + A)} and the normalization of the 
seed photon spectrum Ncomptb (see Appendix A for the definition of Ncomptb) ■ 

In Figure [2] we illustrate our spectral model as a basic model for fitting the BeppoSAX 
and RXTE spectral data for 4U 1728-34. We assume that accretion onto a neutron star 
takes place when the material passing through the main two regions, a geometr ically thin 
accretion disk [standard Shakura-Sunyaev disk, see IShakura fe Sunyaev I (119731 )] and the 
transition layer (TL), where NS and disk soft photons are upscattered off hot electrons. In 
other words, in our picture, the emergent thermal Comptonization spectrum is formed in the 
TL region, where disk BB-like seed photons and neutron star soft photons are upscattered in 
the relatively hot plasma. Some fraction of these seed soft photons can be also seen directly 
by the Earth observer. Red and blue photon trajectories shown in Fig. [2] correspond to soft 
(seed) and hard (up-scattered) photons respectively. 

We show examples of X-ray spectra in Fig. [3] (for BeppoSAX data) and in Figs. HH5] (for 
RXTE data). Spectral analysis of BeppoSAX and RXTE observations indicates that X-ray 
spectra of 4U 1728-34 can be described by the model where its Comptonization component 
can be presented by COMPTB model. Moreover, for broad-band BeppoSAX observations 
this spectra l model component is m odi fied by photoelectric a bsorption at low energies. Also 
following to IDi Salvo et al. I ( 120001 ) and iPiraino et al. I ( 120001 ) suggestions, we add Gaussian 
line at ~6.7 keV and thermal blackbody compone nt at low energies (1 - 4 keV) to improve 



the fit statistics. Note along with these components IDi Salvo et al. I ( 120001 ) included a narrow 
Gaussian line to fit an excess in the residuals around 1.7 keV. We also test the presence of this 
line feature, but the addition of this component to the m odel does not improve a quality of 



the model fit. It is worth noting that ID' Ai et al. I ( 120061 ) analyzed a simultaneous Chandra 



and RXTE observations of the 4U 1728-34 (2002 March, 2-5). They fitted the 1.2 - 35 keV 
continuum spectrum with a blackbody plus a Comptonized component and they fitted large 
residuals at 6 - 10 keV by a broad (FWHM~ 2 keV) Gaussian emission line or, alternatively, 
by two absorption edges asso ciated with low ion ized iron and Fe XXV/XXVI. However, in 
the framework of this model, ID' Ai et al. I (120061 ) found no evidence of broad, or narrow Fe 



3 http://heasarc.gsfc. nasa.gov/docs/software/lheasoft/xanadu/xspec/models/comptb. html 
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K lines, between 6 and 7 keV. But using our model wabs*(blackbody+COMPTB+ Gaussian) 
we found an iron line feature during all BeppoSAX and RXTE observations. 

On the top of Figure [3] we demonstrate the best-fit BeppoSAX spectrum using our 
model and in the bottom right panel we show the best-fit spectrum along with A% for the 
model (reduced x 2 =1.16 for 445 d.o.f). In particularly we find that an addition of the 
soft blackbody-like component of temperature Tbb =0.5—0.7 keV to the model significantly 
improves the fit quality of the BeppoSAX spectra. The line emission is clearly detected in 
the range from 5 to 8 keV as one can see from the left bottom panel of Fig. [3J We show that 
this line is quite broad and it is much wider than the instrumental response whose width is 
smaller than 0.02 keV0. Thus we include a simple Gaussian component whose parameters 
are a centroid line energy Eu ne , the width of the line o~u ne and the normalization Nu ne in the 
model to fit the data in the 6-8 keV range. We also include the interstellar absorption with 
a column density Nh in the model. It should be mentioned that we fixed certain parameters 
of the COMPTB component: 7 = 3 (low energy index of the seed photon spectrum) and 
5 = because we neglect an efficiency of the bulk inflow effect vs the thermal Comptonization 
in the case of NS source 4U 1728-34. 

For the BeppoSAX data (see Tables 1, 3) we find that the spectral index a is of 1.03±0.04 
(or the corresponding photon index r = a + 1 is 2.03±0.04). While the temperature of 
the seed photons T s of the COMPTB component changes from 1.2 to 1.3 keV and color 
temperature of the soft Blackbody component Tbb is around 0.6 keV. 

Unfortunately RXTE detectors cannot provide well calibrated spectra below 3 keV while 
the broad energy band of BeppoSAX telescopes allows us to determine the parameters of 
blackbody components at soft energies. Thus, in order to fit the RXTE data we have to fix 
the temperature of blackbody component at a value of T B b =0.7 keV obtained as an upper 
limit in our analyze of the BeppoSAX data. The best-fit spectral parameters using RXTE 
observations are presented in Table 4. In particular, we find that electron temperature T e of 
the COMPTB component varies from 2.5 to 15 keV, while photon index V is almost constant 
(r = 1.99 ± 0.02 ) for all observations. It is worth noting that the width ou ne of Gaussian 
component does not vary much and it is in the range of 0.3 - 0.6 keV. Color temperature 
T s of COMPTB component is around 1.3 keV which is consiste nt with that using the 



Bep poSAX data set of ou r analysis (Table 3) and previous studies by lDi Salvo et al. I (120001 ) 



and lPiraino et al. I (120001 1 . 



We fixed the value of the COMPTB parameter \og(A) to 1 when the best-fit values of 
\og(A) ^> 1 because in any case of log(A) » 1 a Comptonization fraction / = A/(l + A) 



See ftp://heasarc.gsfc.nasa.gov/sax/cal/responses/98_ll 
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is approximately 1 and variations of A ^> 1 do not improve fit qualit y any more. We use a 
value of hydrogen column Nh = 2.73 x 10 22 cm" 2 , which was found by lPiraino et al. I ( 120001 ) . 
Systematic error of 0.5% has been applied to all analyzed RXTE spectra. 

In Figure H] we show an example of the best-fit RXTE spectrum of 4U 1728-34 for the low 
luminosity state in units of E *F(E) (top) and the spectrum in counts units (bottom panels) 
with Ax for the 30042-03-01-00 observation. On the left bottom panel we demonstrate a 
fit of the model wahs * COMPTB (x 2 ed =2.1 for 61 d.o.f) and on the right, the same as 
the latter one but we add an iron Gaussian line and the blackbody component using the 
model wabs * (blackbody + COMPTB + Gaussian) for which we obtain x 2 ed =1.18 for 57 
d.o.f. The best-fit model parameters for this observation are T=1.99±0.02, T e =10.4±0.3 
keV and -Ez ine =6.54±0.03 keV (see more details in Table 4). Red, violet and blue lines stand 
for blackbody, COMPTB and Gaussian components, respectively. We also apply the same 
procedure to the spectrum during the high luminosity state and in the Figure [5] we present 
the results for the 50023-01-12-00 observation. On the left bottom panel we show a fit of the 
model wabs*COMPTB (% 2 ed =1.79 for 61 dof) and on the right we present the same as the 
latter one but adding an iron Gaussian line and the blackbody components using the model 
wabs * (blackbody + COMPTB + Gaussian) for which x 2 ed =1.2 for 57 d.o.f. The best-fit 
model parameters in this case are r=1.99±0.02, T e =5.5±0.1 keV and _E/j ne =6.75±0.04 keV 
(see more details in Table 4). 

The adopted spectral model shows a very good performance throughout all data sets 
used in our analysis. Namely, a value of reduced Xred = X 2 /^dof, where N do f is a number of 
degree of freedom, is less or around 1.0 for most observations. For a small fraction (less than 
3%) of spectra with high counting statistics Xred reaches 1.5. However, it never exceeds a 
rejection limit of 1.7. Note that the energy range for the cases, in which we obtain the poor 
fit statistic (two among 127 spectra with % 2 =1.7 for 44 dof), are related to iron line region. 
It is possibly that shape of Fe line is more complex than a simple Gaussian (i.e. blend of 
different energies, presence of an edge, or broadening by Comptonization). The fits tend to 
favor a broad line (see Table 4), which might be caused by Comptonization. However, this 
possible complexity is not well constrained by our data. 



M oreover, recent analys is of high-resolution XMM- Newton spectra of 4U 1728-34 Ng et al. 



( 120101 ). lEgron et al. I (120111 )] using different spectral models also reveal evident residuals at 
6-7 keV, which are attributed to the presence of a broad iron emission line. This feature 
can be equally well fitted by a composition of pure iron li ne and co r respo nding; absorp- 
tion edge as well as Laor instead of Gaussian line profile [Ng et al.l ( 120101 )] and also by 
a relativisti cally smeared line c omponent or by a relativistically smeared reflection model 
component Egron et al. I ( )201ll )]. This variety of the line models using in the data analysis 
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for 4U 1728-34 demonstrates complexity of the line appearance in this source. 



It is worth noting that we find some differences between our values of the best-fit 
mo del parameters and tho se in the literature. In particular, the photon index T, estimated 
by iDi Salvo et al. I (120001 ) for observation id= 20674001, is 1.60±0.25 while our value of 
T = 1.9 ±0.2. This difference of the index values can be explained using slightly different 
models. IDi Salvo et al. I (120001 ) included a narrow Gaussian line around 1.7 keV (radiative 
recombination emission from Mg XI) in order to fit an excess of the residuals of the continuum 
model. Howe yer our model r e sult, u sing the BeppoSAX observation (id=20889003), confirms 



the result of iPiraino et al. I ( 120001 ) although we apply a slightly different s pectral models 



Our b est-fit photon index T = 1.9 ± 0.2 are very close to that obtained by IPiraino et al. 



( 120001 ) using the best-fit parameters of COMPTT model [see iTitarchuk I ( 119941 )] electron 
temperature /cT e =3.16±0.03 keV and optical depth (for spherical geometry) r = 11.4 ±0.2. 

Thus using broad band BeppoSAX observations we can well determine all parameters of 
our spectral model while due to the extensive observations of 4U 1728-34 by RXTE we are 
able to investigate the overall pattern of the source behavior during the spectral transitions 
in the 3-60 keV energy range. 



4. Evolution of X-ray spectral properties during spectral state transitions 

We have established common characteristics of the rise-decay spectral transition of 
4U 1728-34 based on their spectral parameter evolution of X-ray emission in the energy 
range from 3 to 60 keV using RXTE/YCAkMEXTE data. In Figures HHS] we present typi- 
cal examples of the RXTE low and high state spectra for 4U 1728-34. In fact, one can clearly 
see from these Figures that the normalization of thermal (blackbody-like) component is a 
factor of 2 higher in the high state than that in the low state, although photon indices T for 
each of these spectra are slightly variable from 1.8 to 2.1 but mostly concentrated around 
T=2 (see that distribution of V on the left-hand panel of Figure [6]). 

We test the hypothesis of r appr 2 using x 2 -statistic criterium. We show the distribu- 
tion of xled(Tappr) = J2iLi {j^^f EEL ^j versus of T appr on the right-hand panel of Figure El 
One can clearly see a sharp minimum of function XredO^appr) around 1 which takes place in 
the range of T appr = 1.99 ± 0.01 with a confidence level 67% and T appr = 1.99 ± 0.02 with a 
confidence level 99% for 127 d.o.f. 

It is important to emphasize that the photon index T is also independent of the luminos- 
ity of blackbody component of COMPTB L^/d\ and the plasma temperature of Compton 
cloud T e (see Figs. EHH}. 
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Using BeppoSXX. data FT11 suggested that the photon index T is about 2 for quite a 
few NS sources which are observed in the different spectral states. FT11 characterize the 
spectral state by a value of electron temperature T e and they show that T = 2 ± 0.2 (or 
a — 1 ± 0.2) when kT e changes from 2.5 to 25 keV. 

A number of X-ray flaring episodes of 4U 1728-34 have been detected with RXTE 
during 1998 - 1999 {R3, R4 sets) and 2000 {R7, R8 sets) with good rise-decay coverage. 
We have searched for common spectral and timing features which can be revealed during 
these spectral transition episodes. We present the combined results of the spectral analysis 
of these observations using our spectral model wabs * {blackbody + COMPTB + Gaussian) 
in Figures [9T-fT0l RXTE/ ASM count rate is shown on the top panel. Further, from the 
top to the bottom, we show the model flux in two energy bands 3—10 keV {blue points) 
and 10 - 60 keV {crimson points). In the next panel we show a change of the transition 
layer electron temperature kT e . One can clearly see the spectral transition from the high 
state to the low state during the time period from MJD 51070 to MJD 51215 when electron 
temperature kT e varies from 3 keV to 15 keV. Normalization of the COMPTB {crimson 
points) and blackbody component {blue points) are shown in the next panel of Figs. [9] and 
[TU1 In particular, one can see from Figures I9HTU1 how COMPTB normalization Nqomptb 
correlates with the variations of ASM count rate and the model flux in 3-10 keV energy band. 
On the other hand, the normalization of the blackbody component Nbb is almost constant 
except at the flaring episode peak, when N BB increases from 0.09 to 0.27 (see blue points in 
Figure [9] at MJD=51093 and 51133). Moreover these flare spectral transitions are related to 
a noticeable increase of soft flux, in the energy range 3-10 keV, and decrease of hard flux, 
that in 10-60 keV (see the second panels from above in Figs. l9l fT0{) . 

The spectral index a {a = T — 1) is presented in the bottom panels of Figures [91— fTTJl 
The index a slightly varies around 1 (or T ~ 2). Photon index T variation over the entire set 
of the observations is presented in Fig. The electron temperature kT e steadily decreases 
during the burst rise (see blue vertical strips in Figs. [9] and ITU]) . Equivalent width of iron 
line EWu ne and the normalizations of COMPTB Ncomptb and blackbody components Nbb 
steadily increase when electron temperature decreases (see Fig. [IT]). 

In fact, a decrease of electron temperature T e of the Compton cloud (transition layer) 
with an increase o f the (disk) soft flux is a well known eff ect and it is explained in details by 



Titarchuk et al. I (119981 ) and iTitarchuk & Fiorito I ( 120041 ) . 



As shown on the right hand panel of Figure [7] the Comptonization fraction / varies from 
0.6 to 0.9. This means that in most cases the soft disk radiation of 4U 1728-34 is subjected to 
reprocessing in Compton cloud and only a small fraction of disk emission component (1 — /) 
is directly seen by the Earth observer. Thus the energy spectrum of 4U 1728-34 during all 
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states is dominated by a Comptonized component while the direct disk emission is always 
weaker and detectable in the flaring episodes only (see Figs. |9 HT0|) . 



Note that for BHs a definition of spectral transition is related to a change of photon 
index T (see e.g. ST09). However there is no one-to-o n e corr espondence between V and cutoff 



or efold) energy Ef id. iTitarchuk fc Shaposhnikov I (120101 ) demonstrate using RXTE data 



for BH binary XTE J1550-564 that Ef id decreases when T increases from 1.4 to 2.1 — 2.2 
until T reaches 2.2 and then increases. Thus for a BH the main parameter used for 

the spectral transition definition is a variable photon index T which monotonically increases 
when a BH source goes to the high state. 

It is important to emphasize once again that in the NS binary 4U 1728-34 the transition 
from the low state to the high state takes place when electron (plasma) temperature changes 
from 15 keV to 2.5 keV, thus following FT11 suggestion, we define the spectral state transi- 
tion in a NS source in terms of the electron temperature T e of the Compton cloud (TL). In 
this case the low state is characterized by high electron temperature T e , while the high state 
is related to low T e . Note that the electron temperature T e is a directly measurable quantity 
and it corresponds to cut-off energy of the spectrum. 

Not all NSs show flares. Only a few NS binaries (such as Z and atoll— sources) display 
spectral transitions during the bursts. Atoll-sources, such as 4U 1728-34, usually show the 
flare transitions. One can establish a substantial difference between a NSs and a BHs due to 
these flare episodes when a source evolves from the low to the high state and when plasma 
temperature of Comptonized region remarkably changes (like in 4U 1728-34 from 2.5 keV to 
15 keV). Thus the NS LMXB 4U 1728-34 shows a steady decrease of electron temperature 
T e from the low state to the high state while the photon index T stays around 2. In contrast, 
in BH sources we see a steady monotonic growth of T which follows by its saturation (see 
e.g. ST09). 



5. Correlations between spectral and timing properties during spectral state 

transitions 

The RXTE light curves have been analyzed using the powspec task from FTOOLS 5.1. 
The timing analysis RXTE/PCA data was performed in 13 - 30 keV energy range using the 
event mode. The time resolution for this mode is 1.2xl0~ 4 s. We generated power density 
spectra (PDS) in 0.1 - 500 Hz frequency range with 0.001-second time resolution. We 
subtracted the contribution due to Poissonian statistics. To model PDS we used QPD/PLT 
plotting package. 
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Previously, timing analysis of 4U 1728-34 was carried out by iDi Salvo et al. I (120011 ) 
as a function of source position in color-color diagram for RXTE data (Rl, Rl sets in 
our terms). In island part of the color-color diagram (corresponding to the hardest energy 
spectra) the power spectrum of 4U 1728-34 shows several features such as a band-limited 
noise component presented up to a few tens of Hz, a low frequency quasi-periodic oscillation 
(LFQPO) at frequences between 20 and 40 Hz, a peaked noise component around 100 Hz 
and one or two QPOs at kHz frequences. During burst evolution (moving along color-color 
diagram) in the lower banana (corresponding to the softer energy spectra) they find a very 
low frequency noise (VLFN) component below ~1 Hz. In upper banana (corresponding to 
the softest energy spectra) the power spectra are dominated by the VLFN with a peaked 
noise component around 20 Hz. 

We find a similar timing behavior of 4U 1728-34 for our data set along with the energy 
spectra. In particular, in Figure [12], we show the details of a typical evolution of X-ray 
timing and spectral characteristics during X-ray flares. The evolution of RXTE/ ASM count 
rate during the 1998 (R3) outburst transition is presented on the top. Red/blue points 
A, B, and C mark moments at MJD = 51122/51128, 51133.27/51133.34 and 51196/51193 
before, during, and after X-ray outburst respectively. In the lower panels (left column) we 
show PDSs for 13-30 keV energy band along with energy spectral diagram E * E(E) (right 
column) for A (top), B (middle) and C (bottom) points of X-ray light curve. The strong 
noise component with a break at 1 - 3 Hz and broad QPOs centered in the range 7-10 
Hz are seen before and after burst but the latter feature is not seen at B moment (see panel 
Bl), i.e. at the X-ray flare peak. During Bl-burst event one can see a noise component with 
the power peak shifted to higher frequency with respect to that at Al/Cl events. In other 
words, the burst power spectrum, in this case, consists of the "high frequency" white-red 
noise component with break at ~40 Hz. On the right hand side we present the E * F(E) 
spectral diagrams (panels A2, B2, C2) related to the corresponding power spectra (panels 
Al, Bl, CI). The data are shown by red points and the spectral model components are 
displayed by blue, black, and dashed purple lines for COMPTB, blackbody and Gaussian 
components respectively. 

Specifically, before the burst (Al red; 30042-03-08-00 RXTE observation, MJD=51122) 
one can see a broken power law noise component with a break at 1 Hz, broad QPO at 
20 Hz (described by Lorentzian with FWHM=11.7±4.5 Hz, rms=6.3±1.0%, x 2 =131 for 
102 d.o.f for 67%confidence level). Later, just before the burst (Al blue; 30042-03-11-00, 
MJD=51128), break frequency of the broken power law noise component shifts from 1 Hz to 
3 Hz and QPO becomes less evident, but still visible at 10—20 Hz range. The low-frequency 
part of v x power diagram below 1 Hz increases right before the burst (Al, blue). During the 
burst (Bl red, 30042-03-14-02, MJD=51133.27; Bl blue, 30042-03-14-01, MJD=51133.34) 
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one can see a white-red noise component with the break frequency shifted to higher frequency 
at about 40 Hz. QPO component is not seen in the power spectrum during the burst at 
frequencies 80 Hz and below. 

We shou ld point out onc e again that the similar behavior of 4U 1728-34 was detected 
previously by lDi Salvo et al. I ( 120011 ) during the 1996—1997 transition from island to banana 
states in col or- color diagram which we call a "burst" transition in this Paper. Note that 
according to iDi Salvo et al. I ( 120011 ) the power spectrum at the upper banana state (at the 
maximum of the burst) consists of two noise components, n amely VLFN (ve r y low frequency 
noise) and HFN (high frequency noise). In addition to IDi Salvo et al. I ( 120011 ) we have 
detected a particular burst when the power spectrum is presented only with a single HFN 
component (see Bl panel of Figure [T2"|) . 

After the outburst (CI blue, 30042-03-18-00, MJD=51193; CI red, 30042-03-20-00, 
MJD=51196) one can see the same features in v x power plot as that before the burst, 
but with slightly different parameters: v\, r ~1, 2 Hz and vqpo ~ 6, 10 Hz (described by 
Lorenzian with FWHM=6.0±2.1/15.0±2.9 Hz, rms=8.1±0.4/10.6±1.4%, \ 2 =139/ US for 
102 d.o.f; all parameter errors correspond to la confidence level) 



In Figure [T3l we present the u x power plot observed on March 7, 2000 (50023-01-01-00 
RXTE observation, MJD=51610) during the quite state in order to compare it with the 
typical sample of PDS during the X-ray flare event (see panel Bl in Fig. [T2"|) . 



6. Discussion 



We show the quasi-constancy of th photon (spectral) index in quite a few observations of 
NS source 4U 1728-34 using BeppoSAX and RXTE observations. In Figures [THE] we present 
plots of the photon index T as functions of the our model parameters: electron temperature 
T e (in keV), the COMPTB normalization (which equals the normalization of NS blackbody 
seed photons) and the Comptonization fraction / = A/(l + A). We obtain these results 
using an application of the first principle spectral model to extensive BeppoSAX and RXTE 
observations of the NS binary source 4U 1728-34. FT11 give an explanation of this index 
stability which was also revealed in other observations of NS binaries. In this Paper, for 
a completeness of our presentation, we review the main points of the FT11 explanation in 
terms of the transition layer model (TLM). 

The energy balance in the transition layer (TL) is dictated by Coulomb collisions 
with protons (gravitational energy release), while inverse Compton and free- free emission 
are the main cooling channels [see a formulation of this problem in the pioneer work by 
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Zel'dovich fc Shakura I (119691 ) and also a similar consideration in iBisnovatyi et al. I (ll980f )]. 



In fact, for the characteristic electron temperature (3 keV < kT e < 30 keV) and density 
values (< 10 -5 g cm -3 ) of these regions in low mass X-ray binaries, Compton cooling domi- 
nates over free-free emission and the relation between the energy flux per unit surface area 
of the coro na Q rm , the radia t ion en ergy density e{r) and electron temperature T e is given 
by [see also iTitarchuk et al. I (119981 )] 



To 



20.2e(r)T e (r) 



(1) 



where tq is Thomson optical depth of the TL. The distribution e(r) is obtained as a solution 
of the diffusion equation 

dh _ 3Q t 
dr 2 ~ 



;tot 



cr 



(2) 



where Q to t = Qcov + Qdisk is the sum of the corona (TL) and intercepted disk fluxes, respec- 
tively. Combination of Eq. (T5]) with two boundary conditions at NS surface and outer TL 
boundary leads us to the formulation of the TL boundary problem (see details in FT11). 
The solution for e(r) is then given by 
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In order to establish the average plasma temperature T e one should estimate the mean energy 
density in the TL as 
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If we now substitute the result of equation (j 
algebra we obtain 

fcT e r (2 + to) 



m P c z 



1 + Qdisk/ Q c 

One should use a formula for spectral index a 



into equation ([T]), after a bit of straightforward 

(5) 



a 



3/9/3 
-2 + Vl + 9' 



(6) 



where 9 = kT e /m e c 2 and /3-parameter defined in ITitarchuk fc Lyubarskijl (119951 ) . If we 
replace (3 by its diffusion limit /Sdiff 



diflf 



r (2 + r ) 



(7) 
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and using equation (jSJ), we obtain the diffusion spectral index as 



3 9 1 + Qdisk/QcOT /o n 

ttdiff = h \ , (o) 

2 V 4 0.25 ' V ' 

or a diff ^1 + 0.8 Qdisk/Qcor and 

T diff « 1 + a dm = 2 + 0.8 g dis k/Qcor (9) 

for Qdisk/ Qcor < 1- 

Thus until Qdisk/Qcor *C 1 photon index r ~ 2. This is precisely that we see in the 
observations of NS 4U 1728-34 (see Figs. M}. 

However in BHs we observe that the photon index monotonically increases with QPO 
frequency and mass ac c retion rate and finally saturates (see ST09, TS09 and ST10). Recently 



Laurent fc Titarchuk I (120111 ). hereafter LT11, have made Monte Carlo simulations of X-ray 
spectral formation in Compton cloud which surrounds a BH and they reproduce the observed 
correlation index vs mass accretion rate. They also demonstrate that the index saturation 
observed in BH sources is a result of two effects, namely cooling of the converging flow by the 
soft disk photons along with the photon capture by a BH. In fact, spectral index is the inverse 
of the Comptonization parameter Y which is proportional to mean number of up-scattering 
N sc and efficiency of up-scattering rj. But in the relatively cold converging flow (CF) photons 
(when mass accretion rate rh in Eddington units is much greater than 1) mostly up-scatter 
off electrons in the direction of the flow for which N sc is proportional to CF optical depth 
tqf (or rh) and rj is inverse proportional to tqf {rh). Thus the spectral (photon) index 
saturates when mass accretion rate increases. This is precisely what was reproduced in the 
MC simulations by LT11. Hence one can conclude that the monotonic growth of the photon 
index T with mass accretion rate followed by its saturation is the observational signature of 
a black hole while the constancy of F (around 2) versus rh (or electron temperature) is the 
NS signature. 



Recently ISoria et. al I (120111 ) found that in ultra luminous X-ray source HLX1 photon 
index changes from 1.8 to 2.95 but they can not find an argument either this source is 
intermediate-mass BH or foreground NS. Comparison of our and FT11 results for NSs and 
ST09 results for BHs (see also LT11) probably indicates that HLX1 is a black hole because 
its photon index changes in the wide range from 1.8 to 2.95 (see Fig. [T4|) but in a NS case 
the index does not vary and has almost constant value around 2 (see Figs. EHH]). 
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7. Summary 

We presented our analysis of the spectral properties observed in X-rays from the neutron 
star X-ray binary 4U 1728-34 during transitions between the low state and the high state. We 
analyze a number of transition episodes from this source observed with Beppo SAX and RXTE 
satellites. For our analysis we use a good spectral coverage and resolution of BeppoSAX 
detectors from 0.1 to 200 keV along with extensive RXTE observations in the energy range 
from 3 to 60 keV. 

We show that the X-ray broad-band energy spectra during all spectral states can be 
adequately fitted by composition of the Blackbody, Comptonized component (COMPTB) 
and Gaussian component. We also show that photon index T of the best-fit Comptonized 
component in 4U 1728-34 is almost constant, about 2 (see Fig. [6]) and consequently almost 
independent of COMPTB normalization L^/D\ Q which is proportional to the disk mass 
accretion rate rh, (see the left hand panel of Fig. [7]) and plasma temperature of Compton 
cloud T e (see Fig. |SJ). Note the soft (disk) photon luminosity L 39 is units 10 39 erg s _1 and 
distance to the source Dio is units of 10 kpc. We should remind a reader that this index 
stability has been recently suggested using a quite a few number of other NS sources, Cyg 
X-2, Sco X-l, GX 17+2, GX 340+0, GX 3+1, GX 349+2, X 1658-298, GS 1826-238, IE 
1724-3045 which were been observed by BeppoSAX at different spectral states [see details in 
FT11]. 

A relatively high value of Comptonized fraction / = 0.6 — 0.9, obtained in the framework 
of our spectral model, indicates to significant reprocessing of X-ray disk emission in Compton 
cloud in 4U 1728-34. 

We also find using BeppoSAX observations that there are two sources of blackbody 
emission, one is presumably related to the accretion disk and another one is related to NS 
surface for which temperatures of soft photons are about 0.7 keV and 1.3 keV, respectively. 

We demonstrate that the photon index T ~ 2 is almost constant when the source 
moves from the low state to the high state, in other words when the plasma temperature of 
Comptonized region varies from 15 to 2.5 keV (see Fig. [8]). 

We present the strong theoretical arguments that the dominance of the energy release 
in the transition layer (TL) with respect to the soft flux coming from the accretion disk, 
Q disk I Q cor <C 1 leads to almost constant photon index r ps 2. 

Thus we argue that this index stability is the intrinsic signature of a NS binary source 
while in BHs the index monotonically changes with mass accretion rate and ultimately sat- 
urates (see ST09). In Figure [T4] we demonstrate the index correlation vs mass accretion 
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rate for a number of BH sources and how the index depends on mass accretion rate in NS 
4U 1728-34. Photon indices of BH candidate sources (GRS 1915+105, GX 339-4, SS 433 
and GRO J1655-40) show clear correlation with mass accretion rate m or with soft photon 
normalization L 3 g / Df Q which is proportional to rh . This correlation is followed by the index 
saturation when rh exceeds a certain level. The behavior of the index for the considered NS 
4U 1728-34 is clear different from that for the sample of BHC sources. 

We acknowledge Chris Shrader and Cristiano Guidorzi for careful reading and editing 
the presented paper. We are very grateful to the referee for his/her valuable comments and 
corrections of the content of this Paper. 



A. On the definition of the normalization of the COMPTB and BMC Models 

The COMPTB and BMC models describe the outgoing spectrum as a convolution of 
the input "seed" blackbody-like spectrum, whose normalizations are Ncomptb an d N BM q 
and color temperature is kT, with the Comptonization Green's function. Similarly to the 
ordinary bbody XSPEC model, bolometrical luminosity 

POO 

Lbol — / ExA(E)dE, (Al) 
Jo 

where A(E) is the photon flux density of blackbody radiation 

E 2 E 

A(E) = 8.052 x K x — — x (exp l)" 1 , (A2) 

and K = Ncomptb N bmc is the normalization of the seed blackbody photon spectrum, 
defined in the same way as the XSPEC bbody model. 

Thus, one can calculate the emergent luminosity of the source as an integral 

f°° z 3 dz 7T 4 
Lboi = 8.052 x K / = 8.052 x K x — . (A3) 



/ e z - 1 15 
On the other hand N BMC , Nqomptb = K = L 39 /Dl Q [s ection 6.2.10 of " User' s Guide of 



X-Ray Spectral Fitting Package XSPEC v.12.6.0" @ and iFarinelli et al. I (|2008)F| , where 



an 

L39 is the source luminosity in units of 10 39 ergs -1 and D\o is the distance to the source 



8 https:/ /heasarc. gsfc.nasa.gov/docs/software/lheasoft/xanadu/xspec/XspecManual.pdf 
9 http:/ /heasarc. gsfc.nasa.gov/docs/xanadu/xspec/models/comptb. html 
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in units of 10 kpc; and if we know K we can find L 3 g having D w . Thus, similarly to the 
ordinary bbody XSPEC model, the normalizations N C omptb and N BMC are a ratio of the 
source (disk) luminosity to the square of the distance 



This implies important property of both COMPTB and BMC models. Namely using 
these models one can correctly evaluate the normalization of the original "seed" component, 
which is presumably a correct mass accretion rate indicator. 
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Table 1. The list of BeppoSAX observations of 4U 1728-34 used in analysis. 



Obs. ID Start time (UT) End time (UT) MJD interval 

20674001 1998 Aug. 23 19:15:27 1998 Aug. 24 09:14:15 51048. 8-51049.4 1 
20889003 1999 Aug. 19 02:01:32 1999 Aug. 20 04:54:32 51409.1-51410.2 2 
Reference flVDi Salvo et al. ~200oll. (2l|Piraino et al. i 2000l 



Table 2. The list of groups of RXTE observation of 4U 1728-34 used in analysis. 



Number of set 


Dates, MJD 


RXTE Proposal ID 


Dates UT 


Rem. 


Ref. 






Rl 


50128-50143 


10073 


Feb. 15 - March 1, 1996 




1, 2, 


3, 


4, 5, 7, 8 


K2 


50710-50728 


20083 


Sept. 19 - Oct. 1, 1997 




3, 4, 


5, 


7, 8 


R3 


51086-51196 


30042 


Sept. 30, 1998 - Jan. 18, 1999 




4 






R4 


51409-51443 


40019 


Aug. 19 - Sept. 22, 1999 


BeppoSAX 


8 






R5 


51237-51359 


40027 


Feb. 27 - Jul. 10, 1999 




6, 8 






R6 


51198-51213 


40033 


Jan. 20 - Feb. 4, 1999 




6, 8 






K7 


51667-51733 


50023 


March 7 - Jul. 8, 2000 




6, 8 






R8 


51652-51657 


50029 


Apr. 18 - 23, 2000 




6, 8 







References: (1) Strohmayer et al. 1996; (2) Ford & van der Klis (1998); (3) van Straaten et al. (2001); (4) Di Salvo et al. 
(2001); (5) Mendez, van der Klis & Ford (2001); (6) Migliari, van der Klis & Fender (2003); (7) Jonker, Mendez & van der 

Klis (2000); (8) TS05 



Table 3. Best-fit parameters of spectral analysis of BeppoSAX observations of 4U 1728-34 in 0.3-60 keV energy 

ranged Parameter errors correspond to la confidence level. 



Observational 


MJD, 


Tbb 






a = 




log(A) 


NcOMPTB 




Np 

line 




X 2 red (d.o.f.) 


ID 


day 


keV 




keV 


r - 1 


keV 






keV 




eV 




20674001 


51048.70 


0.47(3) 


2.65(2) 


1.30(3) 


0.99(7) 


3.76(8) 


0.10(4) 


4.18(3) 


7.4(1) 


0.55(4) 


52(16) 


1.25(457) 


20889003 


51409.50 


0.62(5) 


1.61(1) 


1.21(5) 


1.07(4) 


3.29(4) 


1.06(6) 


3.56(2) 


6.0(1) 


0.43(4) 


51(11) 


1.16(445) 



t The spectral model is wabs * (blackbody + COMPTB + Gaussian), normalization parameters of blackbody and COMPTB components are in units of L37 /d\ 
erg/ s/kpc 2 , where L37 is the source luminosity in units of 10 37 erg/s, dig is the distance to the source in units of 10 kpc and Gaussian component is in units of 

10 -2 X total photons cm~ 2 s~ 1 in line. 



Table 4. Best-fit parameters of spectral analysis of PCA & KEXTE / RXTE observations of 4U 1728-34 in 3-60 keV 
energy ranged The count rate corresponds to 5 PCA units and corrected for background. Parameter errors 

correspond to la confidence level. 



Observational 


Start Time, 


Exposure 


Rate, 


a = 


Te, 


log(A)tt 


Tsrttt 

CO M PT B 


Arttt 

Bbody 


E • 


iv/tt 

(me 


EW|j„ e , 


X 2 a (d.o.f.) 


Fi/Ftttt 


ID 




MJD 




Time, s 


cnt/s 


r- 1 


keV 








keV 




eV 






10073-01 


01-00 


50128 


49 


9632 


1824 


0.98(6) 


3.54(3) 


0.8(1) 


3.4(1) 


2.4(2) 


6.6(1) 


1.2(1) 


100(30) 


1.16(64) 


2.82/1.43 


10073-01 


01-000 


50128 


494 


12180 


1752 


1.00(3) 


3.71(4) 


0.7(1) 


3.4(2) 


2.09(4) 


6.68(8) 


0.7(1) 


106(20) 


1.1(64) 


2.69/1.49 


10073-01 


02-00 


50129 


27 


9488 


2071 


1.00(2) 


2.91(2) 


l.oott 


4.05(6) 


2.85(5) 


6.63(1) 


0.57(1) 


120(20) 


0.98(65) 


3.26/1.31 


10073-01 


-02-000 


50129 


27 


9840 


1927 


1.11(7) 


3.18(7) 


0.94(6) 


3.7(2) 


2.56(2) 


6.64(3) 


0.73(2) 


110(20) 


1.28(64) 


3.03/1.34 


10073-01 


-03-000 


50129 


66 


11550 


2040 


1.02(3) 


3.16(4) 


0.41(3) 


4.53(5) 


3.19(2) 


6.58(2) 


0.32(1) 


150(20) 


1.39(64) 


3.21/1.28 


10073-01 


-03-00 


50129 


93 


9488 


2071 


1.05(4) 


3.17(4) 


0.41(4) 


4.56(7) 


3.28(2) 


6.58(7) 


0.30(1) 


130(10) 


1.15(64) 


3.22/1.27 


10073-01 


-04-000 


50131 


46 


12540 


1493 


0.99(2) 


4.64(3) 


0.2(1) 


3.25(2) 


2.67(2) 


6.53(5) 


0.38(1) 


90(10) 


1.32(64) 


2.08/1.46 


10073-01 


-04-00 


50131 


79 


9664 


1387 


0.95(2) 


5.2(1) 


0.18(5) 


3.12(3) 


2.36(3) 


6.50(3) 


0.32(1) 


90(10) 


1.07(64) 


2.08/1.46, 


10073-01 


06-000 


50135 


48 


4160 


994 


0.91(7) 


8.1(3) 


0.27(2) 


2.3(2) 


1.19(3) 


6.50(2) 


0.30(2) 


50(8) 


1.4(64) 


1.45/1.66J 
1.43/1. 71| 


10073-01 


06-00 


50135 


78 


9600 


1021 


0.90(8) 


8.2(2) 


0.30(2) 


2.34(6) 


1.18(2) 


6.51(2) 


0.31(1) 


40(10) 


1.1(64) 


10073-01 


-08-000 


50136 


88 


14560 


1035 


0.90(9) 


8.6(3) 


0.33(1) 


2.32(4) 


1.05(1) 


6.50(3) 


0.32(1) 


40(10) 


1.12(64) 


1.50/1.81 


10073-01 


07-000 


50136 


89 


16000 


972 


0.86(3) 


8.4(1) 


0.29(2) 


2.25(5) 


1.06(3) 


6.48(5) 


0.31(1) 


40(10) 


0.96(64) 


1.41/1.73 


10073-01 


07-00 


50137 


22 


2992 


1058 


0.98(2) 


8.5(4) 


0.39(3) 


2.43(6) 


1.05(2) 


6.45(4) 


0.33(2) 


40(10) 


1.17(64) 


1.54/1.76 


10073-01 


-08-00 


50138 


06 


7888 


1030 


0.89(9) 


8.1(2) 


0.31(2) 


2.39(3) 


1.07(1) 


6.60(2) 


0.32(1) 


40(10) 


1.56(64) 


1.49/1.75 


10073-01 


09-000 


50138 


86 


16510 


1090 


0.91(8) 


8.4(1) 


0.33(1) 


2.51(2) 


1.08(3) 


6.48(3) 


0.34(1) 


50(10) 


1.54(64) 


1.58/1.89 


10073-01 


-09-00 


50139 


20 


7968 


1095 


0.82(9) 


6.9(1) 


0.18(1) 


2.66(3) 


1.64(2) 


6.62(2) 


0.16(1) 


70(10) 


1.05(64) 


1.73/1.74 


10073-01 


10-01 


50143 


00 


12530 


1191 


0.99(2) 


12.2(3) 


0.84(1) 


2.40(4) 


0.92(3) 


6.61(5) 


0.44(5) 


50(20) 


1.09(64) 


1.49/2.46 


20083-01 


01-00 


50710 


24 


9520 


2013 


1.0(2) 


3.01(4) 


0.59(5) 


3.8(3) 


2.75(6) 


6.6(1) 


0.5(2) 


120(20) 


0.86(57) 


4.65/4.3 


20083-01 


01-01 


50710 


51 


4464 


2043 


0.99(7) 


2.96(5) 


0.4(1) 


5.0(3) 


2.01(2) 


6.7(1) 


1.0(1) 


160(30) 


0.79(57) 


3.33/1.2 


20083-01 


01-02 


50711 


31 


4928 


1998 


1.00(5) 


2.89(3) 


0.6(1) 


4.65(4) 


2.6(3) 


6.7(1) 


0.5(1) 


110(20) 


0.8(57) 


3.25/1.17 


20083-01 


01-020 


50712 


24 


11500 


1183 


0.99(7) 


3.02(4) 


0.3(1) 


4.3(3) 


2.1(1) 


6.58(8) 


0.33(8) 


110(10) 


1.15(57) 


3.03/1.15 


20083-01 


02-000 


50712 


35 


3264 


1908 


1.00(5) 


3.06(2) 


0.28(5) 


4.2(1) 


2.3(1) 


6.5(1) 


0.6(3) 


130(10) 


0.81(57) 


3.09/1.14 


20083-01 


02-01 


50712 


65 


12770 


1657 


1.00(8) 


3.02(3) 


0.3(1) 


4.6(1) 


2.1(2) 


6.6(1) 


0.3(1) 


130(8) 


0.89(57) 


3.27/1.20 


20083-01 


03-00 


50714 


98 


3520 


2281 


0.99(9) 


2.96(2) 


0.5(2) 


5.45(3) 


2.8(6) 


6.7(1) 


0.4(1) 


140(10) 


0.72(57) 


3.71/1.36 


20083-01 


-03-000 


50715 


01 


15650 


2431 


1.12(9) 


2.85(2) 


l.oott 


5.69(8) 


2.19(1) 


6.7(1) 


0.3(1) 


130(10) 


0.93(58) 


3.96/1.44 


20083-01 


04-00 


50717 


51 


16780 


1644 


0.99(4) 


3.22(7) 


0.15(5) 


3.5(3) 


2.2(1) 


6.59(9) 


0.33(9) 


100(10) 


0.92(57) 


2.65/1.01 


20083-01 


04-01 


50718 


38 


16420 


1388 


0.99(2) 


4.21(8) 


0.09(3) 


3.1(1) 


1.85(9) 


6.5(1) 


0.5(3) 


90(10) 


1.39(57) 


2.19/1.10 



Table 4 — Continued 



Observational 


otart 1 lme, 


Exposure 


Rate, 


ct = 


I e , 


1 „„/ A \f t 

log(AJ 


ISjttt 

"COMPTB 


Afttt 

Bbody 




line 




*red (O-O.t.) 


P /ptttt 


ID 


ft J TT~1 

MJU 


Time, s 


cnt /s 


1 1 


kev 








kev 




ev 






20083-01-03-02 


50721.18 


6736 


1002 


0.99(3) 


3.02(3) 


0.28(4) 


3.54(3) 


2.15(6) 




0.4(1) 


100(10) 


0.87(57) 


2.57/0.95 


20083-01-03-020 


50721.19 


12580 


999 


0.99(7) 


2.96(4) 


0.4(1) 


3.51(4) 


2.17(3) 


6.6d) 


0.28(9) 


100(10) 


0.96(57) 


2.56/0.95 


20083-01-04-02 


50722.25 


3360 


936 


0.99(6) 


2.95(4) 


0.8(2) 


3.1(1) 


2.4(3) 


6.7(1) 


0.3fl) 


90(10) 


0.73(57) 


2.39/0.88 


20083-01-04-020 


50722.26 


16540 


1511 


0.99(6) 


3.08(5) 


0.3(1) 


3.4(2) 


2.1(1) 


6.6(1) 


0.3(1) 


100(10) 


0.98(57) 


2.45/0.92 


30042-03-01-00 


51086.30 


8432 


1795 


0.99(2) 


10.4(3) 


0.71(3) 


3.92(5) 


0.99(8) 


6.54(3) 


0.61(1) 


50(20) 


1.21(57) 


2.62/3.67 


30042-03-02-00 


51093.21 


9712 


2421 


0.97(5) 


3.16(4) 


0.41(5) 


5.35(8) 


3.34(4) 


6.53(2) 


0.55(2) 


140(8) 


1.29(57) 


3.90/1.60 


30042-03-01-01 


51109.94 


1920 


1528 


0.94(3) 


11.1(6) 


1.0(2) 


3.19(2) 


1.05(9) 


6.58(8) 


0.53(7) 


30(20) 


1.1(57) 


2.19/3.44 


30042-03-02-01 


51110.01 


2352 


1525 


0.85(9) 


9.6(2) 


0.52(4) 


3.52(3) 


0.93(2) 


6.53(4) 


0.49(3) 


40(20) 


0.96(57) 


2.19/3.39 


30042-03-03-01 


51110.08 


4624 


1584 


0.89(9) 


9.8(2) 


0.60(3) 


3.64(4) 


0.85(2) 


6.55(3) 


0.49(2) 


40(20) 


1.06(57) 


2.28/3.48 


30042-03-04-00 


51113.01 


6608 


1583 


0.85(8) 


9 6(1) 


0.50(2) 


3.85(2) 


1.12(3) 


6.55(2) 


0.49(2) 


40(20) 


1.07(57) 


2.27/3.54 


30042-03-01-03 


51112.94 


2048 


1589 


0.86(9) 


9 7(2) 


0.56(4) 


3.79(5) 


0.89(2) 


6.54(2) 


0.52(3) 


60(20) 


1.22(57) 


2.29/3.516 


30042-03-01-04 


51115.07 


3312 


1599 


0.85(9) 


9.9(2) 


0.53(3) 


3.68(2) 


0.94(3) 


6.53(2) 


0.51(2) 


50(20) 


1.07(57) 


2 29/3 6i 


30042-03-05-00 


51115.14 


6672 


1595 


0.88(7) 


10.2(1) 


0.69(3) 


3.52(4) 


0.92(1) 


6.57(3) 


0.54(2) 


70(20) 


0.93(57) 


2.28/3.610 


30042-03-07-01 


51119.94 


3376 


1719 


0.89(9) 


9.6(2) 


0.53(2) 


4.03(3) 


0.9(2) 


6.50(1) 


0.47(2) 


70(10) 


0.92(57) 


2.48/3.73 


30042-03-07-00 


51120.00 


6752 


1670 


0.88(6) 


9 6(1) 


0.59(2) 


3.92(2) 


0.89(4) 


6.53(4) 


0.51(2) 


70(20) 


1.15(57) 


2.40/3.66 


30042-03-08-00 


51122.17 


5296 


1637 


0.89(7) 


10.10(2) 


0.67(3) 


3.67(3) 


0.92(3) 


6.57(1) 


0.55(2) 


50(30) 


1.09(57) 


2.40/3.66 


30042-03-10-00 


51127.77 


4432 


1787 


0.93(3) 


8.2(2) 


0.46(2) 


4.35(3) 


1.19(3) 


6.51(2) 


0.52(2) 


60(20) 


1.3(57) 


2.62/3.34 


30042-03-10-01 


51 1 27 88 


2656 


1730 


0.92(2) 


7 9C2~I 


0.51(6) 


3 96('8'1 




6.53(4) 


69("9~) 


SOI^O - ) 


88(57~l 


2 54/3 1 7 


30042-03-11-00 


51127.93 


1013 


1752 


0.96(3) 


8.4(3) 


0.55(3) 


4.07(3) 


1.15(2) 


6.46(2) 


0.75(5) 


70(30) 


0.95(57) 


2.58/3.23 


30042-03-12-00 


51128.60 


16640 


1776 


0.98(2) 


7.95(9) 


0.52(2) 


4.14(4) 


1.21(2) 


6.49(4) 


0.69(4) 


60(30) 


0.97(57) 


2.63/3.08 


30042-03-13-00 


51128.93 


9824 


1765 


0.99(1) 


8.13(9) 


0.53(4) 


4.09(2) 


1.23(2) 


6.47(3) 


0.72(5) 


80(30) 


0.76(57) 


2.61/3.08 


30042-03-14-02 


51133.27 


880 


2059 


0.9(1) 


4.4(1) 


0.3(3) 


4.73(3) 


2.71(2) 


6.49(5) 


0.57(7) 


160(20) 


0.93(57) 


3.23/1.88 


30042-03-14-01 


51133.34 


2432 


2033 


1.00(9) 


4.46(7) 


0.3(3) 


4.61(5) 


2.72(9) 


6.62(3) 


0.54(4) 


150(20) 


1.12(57) 


3.18/1.04 


30042-03-14-00 


51133.41 


2160 


2194 


0.99(6) 


4.09(2) 


0.3(2) 


5.12(4) 


2.97(8) 


6.60(3) 


0.62(4) 


160(20) 


1.18(57) 


3.46/1.90 


30042-03-15-00 


51133.55 


13700 


1989 


1.00(3) 


4.79(3) 


0.3(1) 


4.5(2) 


2.71(4) 


6.55(6) 


0.55(2) 


120(20) 


1.3(57) 


3.10/1.90 


30042-03-16-00 


51134.00 


14940 


1787 


1.00(5) 


5.48(3) 


0.35(5) 


4.28(7) 


2.32(4) 


6.54(1) 


0.55(2) 


120(20) 


0.96(57) 


3.75/2.10 


30042-03-17-00 


51134.54 


3200 


1686 


1.00(3) 


6.9(2) 


0.40(8) 


3.91(6) 


1.79(7) 


6.54(4) 


0.55(4) 


120(20) 


1.2(57) 


3.75/2.10 



Table 4 — Continued 



Observational 


otart lime, 


Exposure 


Ftatc, 


OL — 




log(AJ 


lvjttt 

iy COMPTB 


Afttt 

Bbody 


t-'line ) 


wttt 




Xred (d-O.t.) 


V /ptttt 


ID 


1\ ^ TT"\ 
1V1JD 


Time, s 


cnt/s 


1—1 


kev 








kev 




ev 






30042-03-18-00 


51193.32 


9488 


1232 


0.997(9) 


14.4(3) 


1.07(1) 


2.58(3) 


0.97(6) 


6.49(2) 


0.60(3) 


30(20) 


1.12(57) 


1.78/2.79 


30042-03-19-01 


51195.20 


2432 


1253 


1.000(8) 


14.3(5) 


1.01(1) 


2.47(4) 


0.95(9) 


6.53(2) 


0.54(2) 


40(20) 


1.15(57) 


1.82/2.83 


30042-03-19-00 


51195.26 


9568 


1286 


0.99(3) 


14.1(3) 


0.96(5) 


2.58(3) 


0.9(1) 


6.54(1) 


0.55(1) 


50(20) 


0.94(57) 


1.86/2.90 


30042-03-20-00 


51196.94 


4608 


1365 


0.99(4) 


12.06(4) 


0.83(6) 


3.13(2) 


0.90(6) 


6.51(2) 


0.56(2) 


50(10) 


1.21(57) 


1.98/2.95 


40033-06-01-00 


51198.12 


10030 


1124 


0.99(7) 


13.3(3) 


1.00++ 


2.87(3) 


0.86(2) 


6.67(1) 


0.30(4) 


70(8) 


1.29(58) 


1.99/3.20 


40033-06-02-00 


51200.19 


9792 


1219 


0.99(5) 


9.8(2) 


0.61(3) 


3.35(4) 


0.75(4) 


6.63(1) 


0.43(4) 


60(10) 


1.16(57) 


2.20/2.94 


40033-06-02-01 


51201.91 


16420 


1293 


0.99(6) 


7.93(8) 


0.44(2) 


3.89(2) 


1.04(3) 


6.64(2) 


0.40(4) 


60(8) 


1.04(57) 


2.37/2.68 


40033-06-02-03 


51203.98 


13230 


1084 


1.00(3) 


5.44(3) 


0.35(6) 


4.33(8) 


1.98(4) 


6.66(3) 


0.36(1) 


40(10) 


1.07(57) 


2.64/2.08 


40033-06-02-04 


51205.92 


6576 


1685 


0.99(4) 


5.43(5) 


0.33(7) 


4.3(1) 


1.96(5) 


6.65(2) 


0.38(2) 


40(10) 


0.96(57) 


2.63/2.06 


40033-06-02-07 


51206.06 


2784 


1634 


0.99(2) 


5.81(9) 


0.32(9) 


4 2(1) 


1.88(3) 


6.61(3) 


0.35(3) 


40(10) 


0.98(57) 


2.54/2.07 


40033-06-02-05 


51206.12 


2704 


1716 


1.00(3) 


5.37(6) 


4(1) 


4 3(2) 


1.73(8) 


6.64(3) 


0.44(3) 


50(10) 


0.87(57) 


2.68/2.016 


40033-06-02-06 


51206.19 


3120 


996.9 


1.00(4) 


5.55(8) 


4(1) 


3 8(1) 


1.89(3) 


6.69(5) 


0.49(3) 


50(10) 


1.19(57) 


2. 45/1. 


40033-06-03-00 


51207.18 


3264 


1491 


1.00(2) 


6.6(2) 


0.35(7) 


3.69(6) 


1.48(6) 


6.61(3) 


0.39(3) 


30(10) 


1.19(57) 


2.30/2.1|2 


40033-06-03-06 


51207.27 


864 


911.7 


1.00(3) 


6.8(3) 


8(1) 


3 5(1) 


1.33(9) 


6.58(6) 


0.39(6) 


40(10) 


1.2(57) 


2.22/2.10 


40033-06-03-01 


51208.92 


6608 


1207 


0.99(2) 


5.91(7) 


0.31(6) 


3.68(8) 


1 74(6) 


6.69(2) 


0.33(2) 


30(10) 


0.97(57) 


2.28/1.88 


40033-06-03-02 


51210.12 


1856 


1635 


0.99(3) 


4.27(3) 


0.24(6) 


3.81(5) 


2.16(6) 


6.68(3) 


0.41(4) 


40(10) 


1.16(57) 


2.62/1.49 


40033-06-03-07 


51211.78 


2368 


1379 


1.00(2) 


5.46(9) 


0.34(9) 


3.4(1) 


1.98(4) 


6.74(3) 


0.35(3) 


30(10) 


0.93(57) 


2.15/1.64 


40033- 06-03-OS 


51211.91 


12350 


1112 


1.00(1) 


5.71(5) 


30(6) 


3.37(7) 


1.54(3) 


6.68(1) 


0.36(1) 


30(10) 


1.37(57) 


2.10/1.64 


40033-06-03-05 


51213.91 


5152 


1601 


0.99(3) 


3.98(2) 


0.25(6) 


3.74(5) 


2.07(4) 


6.68(2) 


0.39(2) 


30(10) 


1.32(57) 


2.58/1.59 


40027-06-01-01 


51237.04 


4528 


1023 


1.06(2) 


7.8(3) 


0.59(2) 


2.38(4) 


0.78(2) 


6.56(4) 


0.54(4) 


70(10) 


1.17(57) 


1.58/1.57 


40027-06-01-02 


51237.18 


3360 


1087 


0.99(3) 


7.3(2) 


0.38(8) 


2.69(6) 


0.89(4) 


6.57(3) 


0.38(3) 


40(10) 


0.74(57) 


1.67/1.66 


40027-06-01-04 


51238.36 


3392 


851.8 


0.99(1) 


7.9(3) 


0.39(7) 


2.42(4) 


0.84(2) 


6.51(3) 


0.39(3) 


40(10) 


1.09(57) 


1.68/1.65 


40027-06-01-05 


51238.49 


2720 


1044 


0.91(9) 


8.2(4) 


0.36(4) 


2.46(5) 


0.86(2) 


6.58(3) 


0.41(3) 


50(10) 


1.72(57) 


1.57/1.85 


40027-06-01-06 


51238.56 


2720 


1110 


0.97(4) 


8.3(4) 


0.37(6) 


2.63(4) 


0.85(4) 


6.60(3) 


0.32(3) 


40(10) 


0.91(57) 


1.68/1.88 


40027-06-01-03 


51238.76 


9616 


1054 


0.999(9) 


9.6(2) 


0.43(3) 


2.45(3) 


0.73(3) 


6.61(1) 


0.31(1) 


40(10) 


0.95(57) 


1.60/1.93 


40027-06-01-07 


51239.90 


3200 


657.9 


0.99(2) 


8.9(3) 


0.41(7) 


2.41(3) 


0.85(6) 


6.62(3) 


0.30(2) 


40(10) 


1.13(57) 


1.58/1.82 


40027-06-01-08 


51240.04 


1136 


1270 


1.00(4) 


7.1(3) 


0.36(9) 


3.1(1) 


0.94(3) 


6.72(3) 


0.21(3) 


40(10) 


1.09(57) 


1.95/1.99 



Table 4 — Continued 



Observational 


otart lime, 


Exposure 


Rate, 


QL — 


-*e, 


log(A) 


TSjttt 

^COMPTB 


Afttt 

Bbody 




wttt 

line 




X red (d-o.t.) 


V /ptttt 


ID 


1V1JD 


Time, s 


cnt/s 


1—1 


kev 








kev 




e V 






40027-08-01-02 


51359.15 


1056 


934.7 


1.00(4) 


7.2(31 


0.2(1) 


2.89(4) 


1.47(4) 


6.58(5) 


0.29(4) 


40(10) 


0.93(57) 


1.81/1.69 


40019-03-02-02 


51409.17 


1984 


2078 


1.00(3) 


3.02(3) 


0.52(8) 


5.7(2) 


3.16(6) 


6.64(3) 


0.64(3) 


100(20) 


1.37(57) 


4.43/1.71 


40019-03-02-01 


51409.23 


2400 


2102 


0.94(6) 


3.04(6) 


0.32(7) 


6.7(1) 


3.5(3) 


6.67(3) 


0.39(4) 


150(10) 


1.09(57) 


4.47/1.71 


40019-03-02-00 


51409.55 


14370 


1268 


0.94(7) 


2.94(2) 


0.43(5) 


7.54(9) 


3.1(2) 


6.64(2) 


0.59(2) 


170(8) 


1.12(57) 


5.09/1.95 


40019-03-03-01 


51410.16 


2048 


1911 


0.96(8) 


2.87(5) 


0.8(1) 


7.36(7) 


3.2(1) 


6.62(4) 


0.44(2) 


160(10) 


1.07(57) 


5.20/2.18 


40019-03-03-00 


51410.23 


11760 


2321 


0.99(3) 


3.02(2) 


0.40(4) 


7.21(6) 


3.16(4) 


6.62(1) 


0.53(2) 


170(6) 


1.05(57) 


4.95/1.87 


40019-03-01-03 


51442.39 


8800 


438.6 


1.00(6) 


5.7(3) 


0.3(1) 


2.6(1) 


2.7(6) 


6.64(6) 


0.28(2) 


40(10) 


0.83(57) 


1.72/1.35 


40019-03-01-07 


51442.49 


1200 


457.6 


0.99(5) 


5.6(2) 


0.3(1) 


2.6(1) 


2.6(3) 


6.75(6) 


0.22(2) 


40(10) 


1.09(57) 


1.74/1.36 


40019-03-01-04 


51442.60 


1440 


692.9 


0.99(4) 


5.3(2) 


0.3(1) 


2.7(1) 


2.1(4) 


6.73(5) 


0.28(4) 


40(10) 


1.53(57) 


1.79/1.33 


40019-03-01-01 


51442.67 


1152 


705.3 


1.01(6) 


5 2(11 


4(11 


3 0(1) 


2 7(61 


6.69(5) 


0.26(4) 


40(5) 


1.17(57) 


1.83/1.33 


40019-03-01-02 


51442.74 


848 


708.7 


1.00(4) 


5.4(2) 


0.4(2) 


3.0(1) 


5.5(2) 


6.70(7) 


0.21(4) 


40(8) 


1.34(57) 


1.84/1.361 


40019-03-01-00 


51442.80 


12610 


758.2 


0.99(3) 


4.57(8) 


0.26(9) 


3 0(1) 


2.9(2) 


6.73(1) 


0.26(4) 


40(10) 


1.28(57) 


1 99/1 215 


40019-03-01-06 


51443.19 


6288 


759.2 


0.99(2) 


4.80(5) 


0.24(9) 


3.1(1) 


3.4(4) 


6.72(3) 


0.23(2) 


40(10) 


0.91(57) 


2.06/1.331 


50023-01-01-00 


51610.70 


2656 


551.7 


0.95(3) 


10.07(6) 


0.63(6) 


4.04(5) 


0.93(5) 


6.56(3) 


0.46(3) 


50(8) 


1.03(57) 


1.99/2.83 


50023-01-02-00 


51613.69 


2544 


801.5 


0.96(4) 


10.4(2) 


1.00(1) 


2.93(4) 


0.85(4) 


6.62(3) 


0.45(3) 


50(10) 


1.52(57) 


1.98/2.74 


50023-01-03-00 


51615.89 


3376 


788.8 


0.89(6) 


8.5(2) 


0.40(2) 


3.21(3) 


0.9(3) 


6.56(3) 


0.36(3) 


40(81 


1.75(57) 


1.96/2.55 


50023-01-04-00 


51619.81 


3296 


520 


0.93(3) 


7.8(3) 


0.42(4) 


2.96(4) 


1.07(5) 


6.58(3) 


0.35(3) 


40(10) 


1.07(57) 


1.91/2.21 


50023-01-05-00 


51622.81 


3328 


773.5 


0.94(6) 


6 7(31 


0.22(2) 


3.49(3) 


3.69(4) 


5.99(7) 


0.07(1) 


30(10) 


76(571 


1 98/1 89 


50023-01-06-00 


51625.73 


2992 


536.6 


0.98(7) 


5.5(2) 


0.32(5) 


3.11(5) 


1.5(2) 


6.58(4) 


0.31(3) 


40(8) 


1.04(57) 


2.04/1.56 


50023-01-07-00 


51628.85 


2928 


1896 


0.97(9) 


4.1(2) 


0.05(5) 


6.5(1) 


3.9(2) 


6.73(2) 


0.72(3) 


150(6) 


0.88(57) 


4.01/2.02 


50029-23-01-00 


51652.4 


2256 


280.5 


1.00(3) 


13.6(8) 


0.7(1) 


1.4(2) 


0.63(4) 


6.50(5) 


0.29(3) 


40(10) 


1.2(57) 


1.03/1.44 


50029-23-01-01 


51652.46 


2224 


417.3 


0.89(6) 


9.9(5) 


0.34(5) 


1.68(4) 


0.62(6) 


6.54(5) 


0.19(3) 


30(10) 


0.91(57) 


1.04/1.42 


50029-23-02-00 


51657.13 


2464 


906.8 


0.89(3) 


8.7(4) 


0.39(3) 


2.17(5) 


0.68(5) 


6.57(3) 


0.27(3) 


40(10) 


0.95(57) 


1.38/1.77 


50029-23-02-01 


51657.19 


2480 


1027 


0.97(6) 


7.8(3) 


0.46(4) 


2.51(6) 


0.9(1) 


6.63(3) 


0.26(3) 


40(8) 


1.29(57) 


1.57/1.80 


50029-23-02-02 


51657.58 


5664 


773.8 


0.96(2) 


7.9(2) 


0.38(2) 


2.49(3) 


1.2(2) 


6.62(2) 


0.27(2) 


40(10) 


0.94(57) 


1.52/1.72 


50029-23-02-03 


51657.73 


7408 


384.3 


0.99(8) 


8.6(2) 


0.40(4) 


2.28(4) 


1.3(3) 


6.58(2) 


0.26(2) 


40(10) 


0.81(57) 


1.43/1.64 


50023-01-12-00 


51663.96 


1792 


1195 


0.99(3) 


5.5(1) 


0.3(1) 


2.9(1) 


2.6(4) 


6.75(4) 


0.28(3) 


40(10) 


1.2(57) 


1.90/1.43 
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Table 4 — Continued 



Observational 


otart lime, 


Exposure 


Rate, 


Oi — 


T 1 

-*e, 


log(A) 1 i 


Hjttt 

^COMPTB 


jyttt 

Bbody 


inline i 


line 






p, /ptttt 


1JJ 


i\if Tr\ 
MJD 


Time, s 


cnt / s 


1—1 


ke V 








KC V 




e V 






50023-01-13-00 


51667.55 


2592 


1715 


1.02(4") 


2.95(4) 


0.48(9) 


3.8(1) 


2.3(3) 


6.64(3) 


0.36(3) 


50(8) 


1.1(57) 


2.88/1.08 


50023-01-14-00 


51669.55 


2528 


961.9 


1.0(2) 


4.07(9) 


0.2(4) 


3.07(6) 


2.95(9) 


6.72(5) 


0.27(3) 


40(6) 


1.24(57) 


2.01/1.06 


50023-01-15-00 


51673.61 


2048 


313.9 


1.00(4) 


10.1(6) 


0.47(9) 


1.73(3) 


0.74(4) 


6.51(4) 


0.25(4) 


40(10) 


1.09(57) 


1.16/1.42 


50023-01-17-00 


51679.72 


2688 


368.7 


1.01(6) 


9.4(4) 


0.6(1) 


2.15(4) 


0.40(5) 


6.55(4) 


0.32(3) 


50(10) 


0.78(57) 


1.41/1.74 


50023-01-18-00 


51682.51 


2656 


372.2 


0.89(6) 


7 9(3) 


0.32(3) 


2.61(3) 


1.5(4) 


6.63(5) 


0.23(4) 


40(10) 


0.92(57) 


1.55/1.86 


50023-01-19-00 


51685.50 


2720 


632.4 


0.97(5) 


7 3(2) 


0.35(3) 


2.8(3) 


1.9(2) 


6.70(4) 


0.26(3) 


40(6) 


0.81(57) 


1.72/1.76 


50023-01-20-00 


51688.76 


2608 


613.7 


0.94(4) 


8 1(3) 


0.38(3) 


2.67(5) 


1 0(3) 


6.57(4) 


0.28(3) 


50(8) 


0.92(57) 


1.63/1.92 


50023-01-21-00 


51691.71 


3056 


654.4 


1.01(9) 


9.4(4) 


0.62(7) 


2.80(3) 


0.58(7) 


6.56(3) 


0.32(3) 


60(10) 


1.19(57) 


1.73/2.19 


50023-01-22-00 


51695.33 


2720 


476.1 


1.01(9) 


8 6(2) 


0.56(6) 


3.02(4) 


0.81(3) 


6.59(4) 


0.31(4) 


60(10) 


1.1(57) 


1.82/2.22 


50023-01-23-00 


51697.45 


2112 


248 


0.99(2) 


7.8(3) 


0.52(8) 


2.89(6) 


0.92(4) 


6.68(6) 


0.23(5) 


40(10) 


0.89(57) 


1.83/2.14 


50023-01-24-00 


51709.69 


2560 


188.5 


0.97(4) 


9.5(6) 


0.7(2) 


1.96(3) 


0.5(2) 


6.64(5) 


0.30(4) 


40(10) 


1.29(57) 


1.38/1.83 


50023-01-25-00 


51712.35 


2672 


803.3 


0.90(4) 


7.8(2) 


0.30(1) 


2.63(4) 


1.8(5) 


6.72(3) 


0.19(2) 


30(10) 


1.09(57) 


1.56/1.84^ 


50023-01-26-00 


51715.41 


2528 


576.4 


0.99(7) 


13.8(7) 


0.51(4) 


2.36(3) 


0.8(3) 


6.59(4) 


0.16(2) 


30(10) 


1.19(57) 


1.52/2.26 


50023-01-27-00 


51718.93 


3024 


655.6 


0.99(3) 


9.9(5) 


0.7(1) 


2.58(5) 


0.61(6) 


6.67(3) 


0.24(2) 


40(10) 


0.95(57) 


1.72/2.28 


50023-01-28-00 


51721.86 


2880 


238.6 


1.00(2) 


9.3(3) 


0.6(1) 


2.54(4) 


0.61(4) 


6.56(4) 


0.35(2) 


50(6) 


1.07(57) 


1.74/2.23 


50023-01-29-00 


51724.57 


3344 


1063 


0.99(3) 


5.92(9) 


0.27(9) 


3.4(1) 


1.7(2) 


6.57(3) 


0.28(2) 


50(10) 


1.1(57) 


2.13/1.71 


50023-01-30-00 


51727.50 


3088 


950.2 


0.99(2) 


6.4(2) 


0.24(7) 


3.1(1) 


2.1(3) 


6.71(3) 


0.16(2) 


30(10) 


1.18(57) 


1.90/1.62 


50023-01-31-00 


51730.42 


2768 


226 


0.99(4) 


6.7(3) 


0.3(1) 


2.62(8) 


0.9(1) 


6.58(5) 


0.28(5) 


40(10) 


1.02(57) 


1.71/1.59 


50023-01-32-00 


51733.41 


2912 


817.7 


1.0(2) 


3.7(1) 


0.2(1) 


3.4(1) 


1.9(2) 


6.61(4) 


0.35(4) 


50(8) 


1.05(57) 


1.36/1.11 



t The spectral model is wabs * (blackbodv -4- COMPTB + Gaussian), where Nh is fixed at a value 2.73xl0 22 cm" 2 l|Piraino et al. II2000T ): color temperature T« 
and Tbb are fixed at 1.3 and 0.7 keV, respectively (see comments in the text); tt when parameter log(A) 2> 1, it is fixed at 1.0 (see comments in the text), ttt 
normalization parameters of blackbody and COMPTB components are in units of Z/37/cijQ erg/s/kpc 2 , where L37 is the source luminosity in units of 10 37 erg/s, 
dio is the distance to the source in units of 10 kpc and Gaussian component is in units of 10 — 2 X total photons cm~ 2 s~ 1 in line , ^^spectral fluxes (F1/F2) in 

units of xlO -9 ergs/s/cm 2 for (3 - 10) and (10 - 60) keV energy ranges respectively. 
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51000 52000 53000 54000 55000 

MJD, days 

Fig. 1.— Top: Evolution of ASM/RXTE count rate during 1996 - 2006 observations of 
4U 1728-34. Blue rectangles indicate the RXTE/PCA & HEXTE data of pointed observa- 
tions, green triangles show BeppoSAX NFI data used for the analysis. Bottom: ASM/RXTE 
1-day type light curve of 4U 1728-34 during 1996 - 2010. Blue dashed line shows a mean 
count rate and indicates the long-term quasi-periodic variability of mean soft flux during 
~ six years cycle. Green double arrow points out the 1996 - 2000 time interval of RXTE 
observations, used in our analysis (see Rl - R8 intervals). 
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Fig. 2. — A suggested geometry of the system. Disk and neutron star soft photons are 
up-scattered (Comptonized) in the relatively hot plasma of the transition layer (between the 
accretion disk and NS surface). But some fraction of these photons is seen directly by the 
Earth observer. Red and blue photon trajectories correspond to soft and hard (upscattered) 
photons respectively. 
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Fig. 3. — Top : the best-fit spectrum of 4U 1728-34 in E * F(E) units using BeppoSAX 
observation 20889003 carried out on 19th of April 1999. The data are presented by crosses 
and the best-fit spectral model wabs*(blackbody+COMPTB+Gaussian) by green line. The 
model components are shown by red, crimson and blue lines for blackbdody, COMPTB and 
Gaussian components respectively. Bottom panels: the spectrum in units of counts along 
with Ax- Left bottom panel : the best-fit spectrum and Ax for the model fit without the 
line component (reduced x 2= 2.15 for 445 d.o.f) and on right bottom panel: same as that on 
the left one but with an addition of Gaussian (K Q -line) component (reduced x 2= l-16 for 445 
d.o.f). The best-fit model parameters are r=2.07±0.04, T e =3.29±0.04 keV, E line =6.0±0.1 
keV and EWu ne = 51 ± 11 eV (see more details in Table 3). 
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Fig. 4. — The best-fit RXTE spectrum of 4U 1728-34 for the low luminosity state in 
units E * F(E) (top) and the spectra in counts units (bottom panels) with Ax for the 
30042-03-01-00 observation. Left bottom panel: a fit of the model wabs * COMPTB, 
(Xred = 2-1 for 61 d.o.f.) and right bottom panel: the same as the latter one but with ad- 
dition of an iron Gaussian line and the blackbody component, namely using the model 
wabs * (blackbody + COMPTB + Gaussian) (xf ed —l-18 for 57 dof). The best-fit model pa- 
rameters are T=1.99±0.02, T e =10.4±0.3 keV and -E/j ne =6.54±0.03 keV (see more details in 
Table 4). Red, violet and blue lines stand blackbody, COMPTB and Gaussian components 
respectively. 
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Fig. 5. — The best-fit RXTE spectrum of 4U 1728-34 for the high luminosity state 
in units E * F(E) (top) and the spectra in counts units (bottom panel) with Ax for 
the 50023-01-12-00 observation. Left bottom panel: a fit of the model wabs * CompTB 
(xl ed =l-79 for 61 dof) and right bottom panel: the same as the latter one but with ad- 
dition of an iron Gaussian line and the Blackbody component, namely using the model 
wabs * (blackbody + COMPTB + Gaussian) (x 2 red =\fl for 57 dof). The best-fit model 
parameters are r=1.99±0.03, T e =5.5±0.1 keV and -E^ ne =6.75±0.04 keV (see more details 
in Table 4). Red, violet and blue lines stand for Blackbody, COMPTB and Gaussian 
components, respectively. 
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Fig. 6. — Left panel: Histogram (frequency distribution) of the best-fit photon power-law 
index T obtained using a model wabs * (blackbody + COMPTB + Gaussian) for RXTE data 

(1996 - 2000). Right panel: Function x 2 {^a PP r) = ^ Eii ( F< ArT"" ) vs A dashed 

horisontal line indicates the critical residual level Xred = 1 ( see comments in the text). 




0.02 0.04 0.06 0.08 0.2 0.4 0.6 0.8 

COMPTB Normalization, Lg 9 /Df Comptonized fraction, f 



Fig. 7. — On the left photon index T is plotted versus COMPTB normalization Nqomptb = 
L3$/Df and on the right thai versus Comptonized fraction / = A/(l + A) using our spectral 
model wabs * (blackbody + COMPTB + Gaussian) (see details in Tables 3, 4). Blue and 
red points correspond to BeppoSAX and RXTE observations of 4U 1728-34 respectively. 
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Fig. 8. — Photon index plotted versus electron temperature T e (in keV) in the framework 
of our spectral model wabs * (blackbody + COMPTB + Gaussian) (see details in Tables 3, 
4). Blue and red points correspond to BeppoSAX and RXTE observations of 4U 1728-34 
respectively. 
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Fig. 9.— From Top to Bottom: Evolution of the RXTE/ ASM count rate, model flux 
in 3-10 keV and 10-60 keV energy ranges (blue and crimson points respectively), electron 
temperature kT e in keV and COMPTB and blackbody normalizations (crimson and blue 
respectively) during 1998, 1999 flare transition set (R3, R5). In the last bottom panel we 
present an evolution of the spectral index a = V — 1. The rising phases of the burst are 
marked with blue vertical strips. The peak burst times are indicated by the arrows on the 
top of the plot. 
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Fig. 10. — Similar to that presented in Fig. [9] but for the RXTE 2000 flare transition set 
(R7, R8). 
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Fig. 11. — Equivalent width of the iron line (in eV, top), normalizations of COMPTB 
{middle) and Blackbody (bottom) components plotted versus electron temperature T e (in 
keV) in the framework of our spectral model wabs * (blackbody + COMPTB + Gaussian) 
during flare transitions (see also Tables 3, 4). Green points correspond to BeppoSAX obser- 
vations of 4U 1728-34. 
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Fig. 12.- Upper panel: RXTE/ASM count rate during the 1998 - 1999 (R3, R5) out- 
burst transition. Red/blue points A, B, and C mark moments at MJD = 51122/51128, 
51133.27/51133.34 and 51196/51193 (before, during, and after X-ray outburst, respectively). 
Lower panels: PDSs for 13-30 keV energy band (left column) are plotted along with energy 
spectral diagram E * F(E) (right column) related to A , B and C points of X-ray light curve 
(upper panel). The strong noise component related to break at 1 - 3 Hz and broad QPOs 
centered at 7 - 10 Hz are present before and after burst (see panels Al, CI). At the X-ray 
flare peak (see panel Bl) one can see a white-red noise PDS with the break at about 40 Hz. 
On the right hand side panels we present the E * F(E) spectral diagrams (panels A2, B2, 
C2) related to the corresponding power spectra (panels Al, Bl, CI). The data are shown 
by red points and the spectral model components are displayed by blue, black, and dashed 
purple lines for COMPTB, blackbody and Gaussian components respectively. 
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Fig. 13. — The v x power diagram of 4U 1728-34 in 0.1 - 150 Hz range, observed on March 7, 
2000 (50023-01-01-00, MJD=51610). Blue solid line shows the best fit of the power spectrum 
which typically consists of three components: the broad-band noise with break v\> (broken 
power law), low frequ ency QPOs fit by Lorentzians (u s i, i>\) and ~100 Hz bump noise [see 
Pi Salvo et all fboOlh ]. 
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Fig. 14. — Examples of diagrams of photon index E versus mass accretion rate for BH 
candidate sources [GRS 1915+105 ( taken from TS09), GX 339-4 (ST08), SS 433 (ST10) and 
GRO J1655-40 (ST08)] along with that for atoll NS 4U 1728-34. One can see a noticeable 
change of E followed by saturation plateau for BHs as for NS 4U 1728-34 the index slightly 
varies about 2 (see also Fig. |6]). The level for E = 2 is indicated by blue dashed line. 



